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Abstract 

Problem 1.5.7 from Pitman's Saint-Flour lecture notes [9]: Does there exist for 
each n a T-Vj-valucd fragmentation process (II ni fe, 1 < k < n) such that Tl n ,k is distributed 
like the partition generated by cycles of a uniform random permutation of [n] conditioned to 
have k cycles? We show that the answer is yes. We also give a partial extension to general 
exchangeable Gibbs partitions. 



1 Introduction 
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Let [n] = {1, 2, . . . , n}, let P" - be the set of partitions of [n], and let be the set of partitions 
of [n] into precisely k blocks. 

The main result of this paper concerns the partition of [n] induced by the cycles of a uniform 
random permutation of [n]. We begin by putting this in the context of more general Gibbs 
partitions. Suppose that v Ui k, n > 1, 1 < k < n is a triangular array of non- negative reals 
and Wj, j > 1 is a sequence of non-negative reals. For given n, a random partition n is 
said to have the Gibbs(v,w) distribution on V n if, for any 1 < k < n and any partition 
{A 1 ,A 2 , ...,A k }€ VI, we have 



X 
S3 



'{U = {A 1 ,A 2 ,...,A k }) =t;„ 1 jfcJJ«;| j4j 



j'=i 

In order for this to be a well-defined distribution, the weights should satisfy the normalisation 
condition 

n 
k=l 

k 

Bn,*(w) := H w \Aj\- 
{A 1 ,A 2 ,...,A k }eT^j=l 



where 
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B n k(yv) is a partial Bell polynomial in the variables w\, W2, ■ ■ •• Let K n be the number of blocks 
of II. Then it is straightforward to see that 

¥(K n = k) = v njk B n>k (w) 

and so the distribution of II conditioned on the event {K n = k} does not depend on the weights 
v n ,k, n > 1, 1 < k < n: 

F(U = {A 1 ,A 2 ,..., A k }\K n = k) = ll l =1 ™ lA f . 

By a Gibbs(w) partition on V2 we will mean a Gibbs(v, w) partition of [n] conditioned to have 
k blocks. 

For given n, a Gibbs(w) fragmentation process is then a process (H n .k, 1 < k < n) such that 
n ni fc € for all k, which satisfies the following properties: 



(i) for k = 1, 2, . . . , n, we have that H Ut k is a Gibbs(w) partition; 

(ii) for k = 1, 2, . . . , n — 1, the partition II n .fc + i is obtained from the partition n ni fc by splitting 
one of the blocks into two parts. 



If Wj = (j — 1)!, the Gibbs(w) distribution on V k is the distribution of a partition into blocks 
given by the cycles of a uniform random permutation of [n], conditioned to have k cycles. 
Problem 1.5.7 from Pitman's Saint-Flour lecture notes [U] asks whether a Gibbs fragmentation 
exists for these weights; see also [1] for further discussion and for results concerning several 
closely related questions. 

In Section [2] we will show that such a process does indeed exist. Using the Chinese restau- 
rant process construction of a random permutation, we reduce the problem to one concerning 
sequences of independent Bernoulli random variables, conditioned on their sum. In Section 
we describe a more explicit recursive construction of such a fragmentation process. Finally in 
SectionU]we consider the properties of Gibbs(w) partitions with more general weight sequences, 
corresponding to a class of exchangeable partitions of N (and including the two-parameter fam- 
ily of (a, #)-partitions). For these weight sequences we prove that, for fixed n, one can couple 
partitions of [n] conditioned to have k blocks, for 1 < k < n, in such a way that the set of 
elements which are the smallest in their block is increasing in k. This extends the result from 
Section [2] on Bernoulli random variables conditioned on their sum; it is a necessary but not 
sufficient condition for the existence of a fragmentation process. 



2 Existence of a fragmentation process 



Consider the Chinese restaurant process construction of a uniform random permutation of [n], 
due to Dubins and Pitman (see, for example, Pitman [9]). Customers arrive in the restaurant 
one by one. The first customer sits at the first table. Customer i chooses one of the i — 1 
places to the left of an existing customer, or chooses to sit at a new table, each with probability 
So we can represent our uniform random permutation as follows. Let C%, ■ ■ ■ ,C n be a 
sequence of independent random variables such that Cj is uniform on the set 1,2, . . . ,i — 1. 
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Let B\, B2, ■ ■ . , B n be independent Bernoulli random variables, independent of the sequence 
C2,C%,... ,C n and with Bi having mean If Bi = 1 then customer % starts a new table. 
Otherwise, C% gives the label of the customer whom he sits next to. The state of the system 
after n customers have arrived describes a uniform random permutation of [n]; each table in the 
restaurant corresponds to a cycle of the permutation, and the order of customers around the 
table gives the order in the cycle. Write H(Bi, B2, . . . , B n , C2, C3, . . . , C n ) for the random parti- 
tion generated in this way (the blocks of the partition correspond to cycles in the permutation, 
i.e. to tables in the restaurant). 

This construction has two particular features that will be important. Firstly, the number of 
blocks in the partition is simply Y^i=\ S° if we condition on Y17=i = k, we obtain precisely 
the desired distribution on V^, of the partition obtained from a uniform random permutation 
of [n] conditioned to have k cycles. Secondly, if one changes one of the Bi from to 1 (hence 
increasing the sum by 1), this results in one of the blocks of the partition splitting into two 
parts. 

Hence, we can use this representation to construct our sequence of partitions (n n 1 < k < n). 
We will need the following result. 

Proposition 2.1. Let n be fixed. Then there exists a coupling of random variables Bf, 1 < i < 
n, 1 < k < n with the following properties: 



(i) for each k, 



(B 1 ,±S 2 ,..., B n ) — \ ±Si, ±S 2 , . . . , B n 



£^ = fc ; 



i=i 



(ii) for all k and i, if B^ = 1 then B^ +1 = 1, with probability 1. 



So we fix C2, C3, . . . , C n , and define 

n ni fc = H(Bi , B%, . . . , B n , C2, C3, . . . , C n ). 
Then (II n fc, 1 < k < n) is the desired Gibbs fragmentation process. 

It remains to prove Proposition 12.11 Write B k = (B*, B%, . . . , B%), b = (b±, 62, . . . , b n ) and 

p k (b) := P ^B k > n = b) , 

where bi € {0, 1} for 1 < i < n. Clearly, p k (b) is only non-zero for sequences b having exactly k 
l's. Write for the subset of {0, l} n consisting of these sequences and write b -< b' whenever 
b' can be obtained from b by replacing one of the co-ordinates of b which is by a 1. We 
need a process (B 1 , B 2 , . . . , B n ) whose /cth marginal B k has distribution p k and such that, with 
probability 1, B k -< B k+1 for each k. 

It is enough to show that for each k, we can couple B k with distribution p k and B k+l with 
distribution p k+1 in such a way that B k -< B k+1 . (The couplings can then be combined, for 
example in a Markovian way, to give the desired law on the whole sequence) . As noted in Section 
4 of Berestycki and Pitman [1] , a necessary and sufficient condition for such couplings to exist 
with specified marginals and order properties was given by Strassen; see for example Theorem 1 
and Proposition 4 of [8] . (Strassen's theorem may be seen as a version of Hall's marriage theorem 
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Figure 1: The case n = 4 and k = 2, 3. Since Si is always 1, we omit it from the picture. For 
2 < i < 4, a filled circle in position i indicates that = 1; an empty circle indicates that it 
is 0. The arrows indicate which states with k = 3 can be obtained from those with k = 2; the 
numbers are the probabilities Pk(b) of the states. 



[6] and is closely related to the max-flow /min-cut theorem [31 E] ) • The required condition may 
be stated as follows. For C C Sk, write N{C) = {b' € Sk+i : b -<b' some 6 6 C}. We then need 
that for all C C 5^, 

yw< y; (i) 



bee 6'eAr(C) 



Phrasing things a little differently, we need the following proposition. 

Proposition 2.2. Xei Ai, A2, ■ ■ ■ , A m be any collection of distinct k-subsets of [n] and let 

Ej = {B, t = 1 V i G Aj}, l<j<m. 

Then 



F \Eh U E 2 U • • • U E m 
1 < k < n - 1. 



B i = k ) < P ( -Ei U £ 2 U • • • U £ n 



i=l 



i=l 



This is a corollary of the following result from Efron [2] . 

Proposition 2.3. Let <p n : {0, l} n — > ]R+ 6e a function which is increasing in all of its argu- 
ments. Let Li, 1 < i < n, be independent Bernoulli random variables (not necessarily with the 
same parameter). Then 



E 4> n (h,h, ■ ■ ■ ,I n ] 
for allO <k<n - I. 



i=l 



< E 



(f>n(h,h, ■■■ ,In) 



J2L i = k + l 



i=i 



Since Efron's proof does not apply directly to the case of discrete random variables, we give a 
proof here. 



4 



Proof. First let tp : — > R + be a function which is increasing in both arguments. Fix n > 1 
and write X n = £^=1 /i- We will first prove that for < k < n, 

E [i>(X n , I n+1 )\X n + / n+1 = fc] < E [rp{X n , I n +i)\X n + I n+1 =k + l}. (2) 

Let u(i) = P (X n = i) for < i < n. We observe that, as a sum of independent Bernoulli 
random variables, X n is log- concave, that is 

u{i) 2 > u(i - l)u(i + 1), i>l. (3) 

(This follows because any Bernoulli random variable is log-concave, and the sum of independent 
log-concave random variables is itself log-concave, as proved by Hoggar [7]). 

Since ^(0, 0) < ^(1, 0) and ^(0, 0) < ^(0, 1), we have 

K[^(X n ,I n+1 )\X n + I n+1 = 0] =v(o,o) 

p n+ iu(0)^(l, 0) + (1 - p n+ i)n(l)V(0, 1) 



< 



p n+1 u(0) + (1 -p n+ i)u(l) 
= E[i;(X n ,I n+1 )\X n + I n+1 = 1]. 

Similarly, since ip{n — 1, 1) < ?/>(n, 1) and V>( n >0) < ift(n, 1), we have 

p n +iu(n - l)ip(n - 1, 1) + (1 - p n+1 )u(n)ijj(n, 0) 



E [tp(X n , I n+ i)\X n + I n+1 = n] 



p n +iu(n - 1) + (1 - p n+ i)u{n) 

<V>M) 

= E [i/>(X n , I n +i)\X n + I n+1 = n + 1] . 

Suppose now that 1 < k < n — 1. Then 

E [^(X n ,/ n+1 )|X n + J n+1 = fc + 1] - E [V>(X„, J n+ i)|X n + I n+1 = k] 
= p n +Mk)jj(k, !) + (!- p n +i)u(k + l)^(fc + 1, 0) 
p n+1 u(k) + (1 - p n+ i)u(k + 1) 
_ Pn+m(fc - IMA; - 1, 1) + (1 - p n+ i)^(fc)V>(fc, 0) 

Pn+iu(k - 1) + (1 -p n+1 )u(fc) 
= p% +1 u(fc)tt(fc - l)[l{j(k, 1) - ^(fc - 1, 1)] 
d 

(1 - p n+1 ) 2 u{k)u{k + l)[^(fc + 1,0)- i/>(k, 0)] 
d 

p ra+ i(l - p n +i)[u(k - l)u{k + l)V;(fc + 1,0)- u(k) V(fc, 0)] 

p ra+ i(l - p n+ i)[n(fc) 2 ^(fc, 1) - u{k - l)u{k + l)y>(fc - 1, 1)] 
+ d 

where the denominator d is given by 

d = [p n+1 u(k) + (1 - p n+ i)u{k + 1)] [p n+1 u(k - 1) + (1 - p n+ i)u(k)] 

and is clearly non-negative. The first two terms in (jl]) are non-negative because ip is increasing. 
The sum of the third and fourth terms is bounded below by 

Pn+lfl-Ptt+l) 2 _ ^ k _ k + 1 )]^( A , ) x) _ ^( fcj0 )]. 



(4) 
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This is non- negative by the log-concavity property ©. Putting all of this together, we see that 
E ty(X n ,I n+1 )\X n + I n+1 = k]<E [i/>(X n , I n+l )\X n + I n+1 = k + 1] , 

as required. 

We now proceed by induction on n. Note that ([2]) with n = 1 gives the base case. Now define 



'n-l 



4> n ^hJn = E 



<i=l 



f >n(h,h, ■ ■ ■ ,I n ] 



n-l 



j=l 



Assume that we have 



E 



-i(h, h, ■ ■ ■ , In-l) 



n-l 



2> = * 



i=l 



< E 



-i\h,h, ■ ■ ■ , In-i) 



n-l 



i=l 



for < k < n — 2. By this induction hypothesis, ip n \Y^i = i Ii,In) is increasing in its first 

argument. By the assumption that cp n is increasing, ip n (^27=1 In) is increasing in its second 
argument. So by ([2]), 



E 



'n-l 



J=l 



n-l 



h + In = k 



1=1 



< E 



'n-l 



d=l 



n-l 



J2 h + In = k + 1 



i=l 



for < k < n — 1. But by the tower law, this says exactly that 



, Mil, h,-.., In) 

The result follows by induction. 



i=l 



< E 



<j>n(h,h, ■■■ ,In) 



Y J k = k + i 



i=i 



□ 



Proof of Proposition EH Let 4> n (B 1 ,B 2 , ... , B n ) = l{£ 1 u£ 2 u---u£ m }- Then cf> n is increasing in 
all of its arguments and so by the previous proposition we have 



£iU£ 2 U'-'U£„ 



J2B i = k) =E 



i=i 



< E 



t>n(Xl,X2, . . . , X n ) 
1>n(Xi,X2, . . • , X n ) 

E 1 UE 2 U---UE m 



J2B i = k 

i=l 
n 

J2 B i = k +l 



i=l 



J2B i = k + l 



i=l 



□ 



Note that we have only proved the existence of a coupling. There is no reason why it should be 
unique. Indeed, in general, there is a simplex of solutions. For the example given in Figure [H 
the extremes of the one-parameter family of solutions are shown in Figure [2j 
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Figure 2: The two extreme solutions to the coupling problem for n = 4 and k = 2, 3 (as drawn 
in Figure [TJ). The arrows are labelled with the joint probabilities P (B k = b,B k+1 = b'). 

3 A recursive construction 



In this section, we describe a more explicit construction of the Gibbs fragmentation processes, 
which is recursive in n and possesses a certain consistency property as n varies. 

The basic principle is the following simple observation. Suppose we want to create a uniform 
random permutation of [n] conditioned to have k cycles. Then n either forms a singleton, or 
is contained in some cycle with other individuals. If it forms a singleton, then the rest of the 
permutation is a uniform random permutation of [n — 1], conditioned to have k — 1 cycles. If, 
on the other hand, n is not a singleton, then we take a uniform random permutation of [n — 1] 
into k cycles and insert n into a uniformly chosen position. 

So we proceed as follows. A Gibbs fragmentation process on [n] for n = 1 or n = 2 is trivial. 



Suppose we have constructed a process (II"~ ,Il2 _ , . . . ,11™^) on [n — 1], with the required 
marginal distributions and splitting properties. We will derive from it a process (II", IE? , . . . , II™) 
on [n]. 

For each k, the partition II n of [n] into k parts will come from either 



(a) adding a singleton block {n} to II?,; or 

(b) adding the element n to one of the blocks of II^ _1 , by choosing an element C n uniformly 
at random from [n — 1] and putting n in the same block as C n . 



Note that 



' ({n} is a singleton in II n } 



B n = l 



J2Bi = k 



i=l 



which is increasing in k, by the monotonicity results proved in the previous section. So let K n 
be a random variable whose distribution is given by 



1 (K n < k) = P ({n} is a singleton in Il£) ; 
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(1,2,3,4,5,6) {1,2,3,4,5}, {6} 

{1,2,3,4,5} {1,3,5}, {2,4) 

{1,2,3,4} {1,3}, {2,4} 



-i 

(1,2,3} 

~t v'~ 

{1,2} 

~\ 



fl,3},{2} 



(1},{2) 



{1,3,5},(2,4},{6} 
{1,3,5}, {2},{4} 
{1,3},{2},{4} 



[1}.{2},{3} 



{1,3,5},{2},{4},{6} {1,3},{2},{4},{5},{6} 

{1,3},{2},{4},{5} {D,{2},{3},{4},{5} 

—/— 
(1},{2},{3},{4) 



1},{2},{3},{4},{5},{6} 



Figure 3: A possible realisation of the coupling. 

now we take choice (a) above for all k > K n , and choice (b) for all k < K n (using the same 
value of C n for each such k). 

A possible realisation of this recursive process is illustrated in Figure [3l 



4 Partial extension to exchangeable Gibbs partitions 



At the beginning of this paper, we introduced the notion of a partition of [n] with Gibbs(v, w) 
distribution. In general, there is no reason why these partitions should be consistent as n varies. 
That is, it is not necessarily the case that taking a Gibbs(v, w) partition of [n + 1] and deleting 
n + 1 gives rise to a Gibbs(v, w) partition of [n]. If, however, this is the case, we can define an 
exchangeable random partition II of N as the limit of the projections onto [n] as n — > 00. In 
this case, we refer to II as an exchangeable Gibbs(v, w) partition. 

An important subfamily of the exchangeable Gibbs(v, w) partitions are the (a, #)-partitions, 
whose asymptotic frequencies have the Poisson-Dirichlet(a, 6) distribution, for < a < 1 and 
6 > —a, or a < and 6 = m\a\, some m G N (see Pitman and Yor [10] or Pitman [9] for a 
wealth of information about these distributions). Here, the corresponding weight sequences are 



(1 - a)j_i T i, v n>k 



{B + a) { 



fc-l)ta 



+ %- 



i)U 



where (x) m ^p := njliC 3 - + 0' — an d ( x )ow := 1- I n the first P ar t of this paper, we have 
treated the case a = 0: a (0, repartition of [n] has the same distribution as the partition 
derived from the cycles of a uniform random permutation. In view of the fact that the array 
v does not influence the partition conditioned to have k blocks, we have also proved that a 
fragmentation process having the distribution at time A; of a (0, #)-partition conditioned to have 
k blocks exists. The Gibbs(w) distribution corresponding to the case a = — 1 (i.e. to weight 
sequence Wj = j\) is discussed in Berestycki and Pitman pQ; in particular, it is known that a 
Gibbs(w) fragmentation exists in this case. 
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Gnedin and Pitman [5] have proved that the exchangeable Gibbs(v, w) partitions all have w- 
sequences of the form 

Wj = (1 - a)j-xii 

for some — oo < a < 1, where for a = — oo the weight sequence is interpreted as being identically 
equal to 1 for all j. Furthermore, the array v must solve the backward recursion 

V n ,k = ln,kV n +l,k + V n+1>k+1 , 1 < k < n, 

where v\ t \ = 1 and 

{n — ak if — oo < a < 1 
k it a = —oo. 

The case a = 1 corresponds to the trivial partition into singletons and will not be discussed 
any further. 

It seems natural to consider the question of whether Gibbs fragmentations exist for other weight 
sequences falling into the exchangeable class. We will give here a partial extension of our results 
to the case of a general a € [— oo, 1). 

Fix n and consider a Gibbs(v, w) partition of [re] with Wj = (1 — a),_ifi for some a € [— oo, 1). 
Let Bi be the indicator function of the event that i is the smallest element in its block, for 
1 < % < n, and let K n = Y^l=i Bi, the number of blocks. Then we have the following extension 
of the earlier Proposition 12.11 



Proposition 4.1. For each n, there exists a random sequence of vectors ( B±' n , . . . , Bn' n 
k = 1,2, ... ,n such that for each k, 

Bi' n , . . . ,B^ ,n ^j = (Si, . . . ,B n \K n = k) , 
and such that R n ' fe+i > R n ' fc f or a ll k. 



This is precisely what we proved earlier for the case a = (and, indeed, for more general se- 
quences of independent Bernoulli random variables). Namely we show that Gibbs(w) partitions 
of [n], conditioned to have k blocks, can be coupled over k in such a way that the set of elements 
which are the smallest in their block is increasing in k. In the context of random permutations, 
this was enough to prove the existence of a coupling of partitions with the desired fragmentation 
property, using the fact that the random variables Cj in the Chinese restaurant process, which 
govern the table joined by each arriving customer, were independent of each other and from the 
random variables Sj. For this is no longer the case and we cannot deduce the full result 

for partitions from this result for "increments" or "records". In fact, this stronger result is not 
always true, for example when a = — oo and so Wj = 1 for all j. In this case it is known that 
no Gibbs fragmentation process exists for n = 20 and for all large enough n (see for example 
PQ for a discussion). By a continuity argument, one can show similarly that if —a is sufficiently 
large, then for certain n no Gibbs fragmentation process exists. 

As in Section EJ Proposition 14.11 can be proved by induction over n. Suppose we have carried 
out the construction for n — 1, and wish to extend to n. Conditioned on K n = k, we need to 
consider two cases: either B n = 1 and K n _\ = k — 1, or B n = and K n = k. Depending on 
which of these cases we choose, we will set either 

( u n ' k T3n,k\ ( r>n— l,k— 1 D n-l,fe-l -, \ 

[B 1 ,...,B n ' j = \B 1 B n ^ , 1 J 
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or 

B? k ,...,B% k ) = (^-^,...,^1^,0). 

Precisely as in Section O this can be made to work successfully provided that the following 
lemma holds. 

Lemma 4.2. For any n, the probability ¥ (B n = l\K n = k) is increasing in k. 

Equivalently, we are showing that the expected number of singletons in U n , conditioned on 
K n = k, is increasing in k (since, by exchangability, the probability that n is a singleton in IT n 
is the same as the probability that i is a singleton, for any i £ [n]). 

The rest of this section is devoted to the proof of this lemma. 

Define S a (n,k) = Z? nj fc(w). Then 

F(K n = k)=v n<k S a (n,k), (6) 

where the S a obey the recursion 

S a {n + 1, k) = j nt kS a (n, k) + S a (n, k - 1), (7) 

with boundary conditions S a (l,l) = 1, S a (n,0) = for all n, and S a (n,n + 1)=0 for all n. 
These are generalized Stirling numbers. We have 

¥(B n = l\K n = k) = P(J^ n _i = k - l\K n = k) = 

o a {n, k) 

(see Section 3 of [5] for further details). Using the recursion, this is equal to 

i 7n-l,fcffq(n - l,k) \ _1 
S a (n — 1, k — 1) 

For this to be increasing in k, it is equivalent that 

S a (n — 1, k — 1) 



7 n _i >fc 5 Q (n - l,k) 
should be increasing in k. This is implied by the following proposition. 
Lemma 4.3. For all a < 1 and all n and k, 

7n,kS a {n, k) 2 > 7„ > fc + iS' Q ,(ra, k + l)S a (n, k - 1). 

Proof. If < a < 1 then it is sufficient to prove the statement that (S a (n, k), < k < n) is 
log-concave, that is 

S a (n, k) 2 > S a (n, k — l)S a (n, k + 1), 

since in that case 7 n ,fc, defined at ([5]), is decreasing in k. Theorem 1 of Sagan [11] states that 
whenever t n f. is a triangular array satisfying 

tn,k — Cn,fc^n— I dn,ktn— l,k 

for all n > 1, where t n .fc, c n ^ and d n ^ are all integers and such that 
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• c Uy k and d n k are log-concave in k, 

• c ni k-id n ^+i + c n ^+\d n ^-i < ^Cn,kd n> k f° r an n > 1> 

then i raj fc is log-concave in /c. These conditions are clearly satisfied for the generalized Stirling 
numbers S a (n,k), with the exception that the sequence d n ^ is not integer- valued. However, 
Sagan's argument extends immediately to this case also, and so we will not give a proof here. 

We turn now to the case — oo < a < 0. We proceed by induction on n. For n = 2 the statement 
is trivial. Suppose that we have 

(n — ak)S a (n, k) 2 > (n — a(k + l))S a (n, k — l)S a (n, k + 1), 

for < k < n. Using the recurrence 

S a (n + 1, k) = (n — ak)S a (n, k) + S a (n, k — 1), 

we obtain that 

(n + 1 - ak)S a (n + 1, kf - (n + 1 - a(k + l))S a (n + l,k- l)S a (n + 1, k + 1) 
= | (n + 1 — ak){n — ak) 2 S a (n, k) 2 

- (n + 1 - a(k + l))(n - a(/c - l))(n - a(& + l))S' a (n, fe - l)5 a (n, A; + 1) j 

+ |2(n + 1 — ak)(n — ak)S a (n, k)S a (n, k — 1) 

- (n + 1 - a(k + l))(n - a(k + A; + l)S a (n, k-2) 

- (n + 1 — a(/c + l))(n — a(A: — l))S a (n, k)S a (n, k — 1)| 

+ |(n + 1 - ak)S a (n, k - l) 2 - (n + 1 - a(fc + l))5 Q (n, fc - 2)5 a (n, (8) 

We take each of the three terms in braces separately. For the first term, we note that 

(n + 1 — ak)(n — ak) — (n + 1 — a(k + l))(n — a(k — 1)) = a(a — 1). 

Hence, by the induction hypothesis, the first term is greater than or equal to 

a(a — l)(n — ka)S a (n, k) 2 . 

Since a < 0, this is non-negative. For the second term in ([8]), we note that applying the 
induction hypothesis twice entails that 

(n — a(k — l))S a (n, k)S a (n, k — 1) > (n — a(k + l))S a (n, k + l)S a (n, k — 2). 

So 

2(n + 1 - a(k + l))(n - a(k - l))S a (n, k)S a (n, k - 1) 
> (n + 1 - a(k + l))(n - a(k + l))S a (n, k + l)S a (n, k-2) 
+ (n + 1 — a(k + l))(n — a{k — l))S a (n, k)S a (n, k — 1). 

It follows that the second term in braces is bounded below by 

2S a (n, k)S a (n, k — 1) + 1 — ak)(n — ak) — (n + 1 — a(k + l))(n — a(k — 1))] 
= 2a(a — l)S a (n, k)S a (n, k — 1). 
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This is non-negative. Finally, we turn to the third term in braces in ([5]). This is equal to 

[(to — a(k — l))S a (n, k — l) 2 — (n — ak)S a (n, k — 2)S a (n, k)} 
— (a — l)[S a (n, k — l) 2 — S a (n, k — 2)S a (n, k)]. 

By the induction hypothesis (at n and k — 1), the first term is non-negative and so this quantity 
is bounded below by 

(1 — a)[S a (n, k — l) 2 — S a (n, k — 2)S a (n, k)]. 

That S a (n,k — l) 2 — S a {n,k — 2)S a (n,k) > is implied by the induction hypothesis (in fact 
it is a weaker statement) and so, since a < 0, the above quantity is non-negative. So <jSJ) is 
non- negative, as required. The case a = — oo follows similarly. □ 
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